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ABSTRACT 

We use high spatial resolution observations of CO to systematically measure the resolved size-line 
widths luminosity-line width, luminosity- size, and the mass-luminosity relations of Giant Molecular 
Clouds (GMCs) in a variety of extragalactic systems. Although the data are heterogeneous we analyze 
them in a consistent manner to remove the biases introduced by limited sensitivity and resolution, 
thus obtaining reliable sizes, velocity dispersions, and luminosities. We compare the results obtained 
in dwarf galaxies with those from the Local Group spiral galaxies. We find that extragalactic CMC 
properties measured across a wide range of environments are very much compatible with those in 
the Galaxy. The property that shows the largest variability is their resolved brightness temperature, 
although even that is similar to the average Galactic value in most sources. We use these results to 
investigate metallicity trends in the cloud average column density and virial CO-to-H2 factor. We find 
that these measurements do not accord with simple predictions from photoionization-regulated star 
formation theory, although this could be due to the fact that we do not sample small enough spatial 
scales or the full gravitational potential of the molecular cloud. We also find that the virial CO-to-Il2 
conversion factor in CO-bright GMCs is very similar to Galactic, and that the excursions do not show 
a measurable metallicity trend. We contrast these results with estimates of molecular mass based on 
far-infrared measurements obtained for the Small Magellanic Cloud, which systematically yield larger 
masses, and interpret this discrepancy as arising from large H2 envelopes that surround the CO-bright 
cores. We conclude that GMCs identified on the basis of their CO emission are a unique class of object 
that exhibit a remarkably uniform set of properties from galaxy to galaxy. 
Subject headings: galaxies: ISM — ISM: clouds — galaxies: dwarf — galaxies: individual (IC 10, 

M 31, M 33, NGC 185, NGC 205, NGC 1569, NGC 2976, NGC 3077, NGC 4214, 

NGC 4449, NGC 4605) — Magellanic Clouds 



1. INTRODUCTION 

Giant Molecular Clouds (GMCs) are the major reser- 
voirs of molecular gas and the sites of most star forma- 
tion in our Galaxy and other galaxies. Their properties 
set the initial conditions for protostellar collapse, and 
may play a role in determining th e stellar initial mass 
function (Mc Kee fc Ostriken l2007f ) . Moreover, because 
GMCs provide the bulk of the material for forming new 
stars their creation may be the limiting mechanism that 
regulates the rate of star formation in galaxies. There- 
fore, increasing our understanding of their properties and 
distribution throughout the different environments of ex- 
ternal galaxies is likely to provide further insights into 
CMC and stellar formation processes. There is a limited 
amount of information, however, that can be gained from 
studies that resolve the general distribution of molecu- 
lar gas but not the individual molecular clouds. Resolv- 
ing GMCs — to measure their sizes, velocity dispersions, 
and luminosities — is a critical step in understanding the 
processes that ultimately drive galaxy evolution. 



Studies of resolved molecular clouds in the Milky 
Way find that GMCs are in approximate virial equi- 
librium and obe y scaling re lations, commonly known 
as Larson laws (jLarsonl 1198 1). that have their origin 
in the ch aracter of the turbulence in the interstellar 
medi um (|Elmegreen fc Scald l2004 iMcKee fc Ostriked 
|2007| ) . Large scale surveys of the Milky Way show that 
GMCs follow uniform scaling rela tions (jSolomon et al.l 
Il987t lElmegreen fc Falgarondll996l ) , with few differences 
prese nt between those located in the i nner and the outer 
d isk (iHeyer. Carpenter, fc Snel]|[200ll ). 

iLarsonl ( 19811 ) established that velocity dispersion, size, 
and luminosity are correlated in Milky Way GMCs. Ob- 
servations indicate th at GMC line w idths increase as a 
power of their radius (jSolomon et al.lil987) . such that 



a^ « 0.72 i?°-^ km s 
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where cr^ is the one-dimensional velocity dispersion of the 
GMC, and R is its radius measured in parsecs. Equation 
[T]is known as the size-line width relation, and it has been 
shown to hold even within GMCs down to very small 
scales (Heycr fc Brunt 2004; Rosolowsky et al. 2008). 

The size-line width relation is an expression of the 
equilibrium turbulence conditions in the molecular ISM. 
Detailed modeling of GMC line profiles shows that the 
emission is macroturbulent (i.e., the scale size of the 
turbulence is larger than the photon mean free path), 
corresponding to many optically thick clumps that have 
a clump-to-clump velocity dispersion similar to the ob- 
served line width (jWolfire et al.l [19931 ). It is known that 



the observed line widths in all but a few very com- 
pact, quiescent clouds are too large to be thermal; at 
the typical GMC temperatures of 15 — 25 K the ther- 
mal CO velocity dispersion would be < 0.1 km s^^. 
They are rather due to supersonic turbulence within 
the clouds. The dominant sources of turbulence in the 
molecular ISM remain somewhat controversial. For ex- 
ample, it is not yet known whether turbulence in molec- 
ular clouds is primarily intern ally or externally driven 
(e.s... iMcKee k Ostrikeijl2007f) . Contributors must in- 
clude star formation — in the form of jets and winds 
for low mass stars, and winds and expanding HIT re- 
gions for massive stars — and star destruction — in 
the form of supernovae shocks and expanding superbub- 
bles. Global processes that couple the large reservoirs 
of thermal or rotational energy to the turbulent cascade 
— such as spiral shocks and several types of instabil- 
ities that involve shear, magnetic fields, and self grav- 
ity — are important on the largest spatial scales. The 
fact th at turb ulence is present on scales > 100 pc (e.g., 
iBrunti 2003 : Dib fc Burkcrt 2005), and that the dissipa- 
tion lifetime of the largest ed dies is < 10^ yr (|Flecklll98ll 
iStone. Ostriker. fc Gammi"3ll998.) . point to constant (or 
at least frequent) energy injection on large scales. This in 
turn suggests that massive stellar death and gobal mech- 
anisms are predominantly responsible for the bulk of tur- 
bulent energy injection in the Milky Way, altho ugh their 
respective dominance is a matter o f debate (jMacLowl 
I200I IPiontek fc OstrikeJ[2f)0l [2005h . 

The last two Larson relations, the luminosity-line width 
and luminosity-size relations, describe correlations be- 
tween cloud luminosity, Leo, and either velocity disper- 
sion or size (note, as is frequently pointed out, that only 
two of the three Larson re lations are independent). They 
are (jSolomon et al.lll987f) 



and 



ico~ 130 cr^ Kknis-i pc^, 



ico«25i?^-^ Kkms-i pc^, 
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with (7y and R expressed in their corresponding units of 
km s""'^ and pc. 

Thus, the observations show that the virial mass, Mvir, 
relates to the CO luminosity (which is proportional to the 
inferred luminous ma ss Mium) in the following manner 
(|Solomon et al.lll987[) . 
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where ico is expressed in K km s"^ pc^. We will 
call this equation the Galactic mass-luminosity relation. 
Note that the exponent of Leo is some what uncertain , 
and we take here the value reported by ISolomon et al.l . 
Although the relation may not be precisely linear it is 
almost so, and simple model s of molecular clouds can ex- 
plain the observed linearity (jDickman. Snell. fc Schloerbl 
Il986[ ). This approximate proportionality between Leo 
and Mmr within the Galaxy allows the definition of an 
empirical CO-to-H2 conversion factor, XeOi that relates 
the intensity of a GMC in the ^^CO J == 1 -^ transi- 
tion to its H2 column density (A^(H2) = -'i^coleo) and 
ultimately it s mass. For GM Cs around a median mass of 
5 X 10'' M0 ISolomon et al.l obtained a CO-to-H2 factor 



that, when adjusted to the currently accepted distance of 
8.5 kpc to the Galac tic Center, yields Xga » 1-9 x 10^° 
cm-2(K km s'^)"^ (jMcKee fc Ostriked BoOTJ) . In the 
Milky Way, this "virial mass" approach yields conversion 
factors similar to those obtained from non-dynamical 
measurements, strongly suggesting that clouds are in a 
dynamical state intermediate between marginal gravi- 
tational binding and virial equilibrium (e.g., 7-ray re- 
sults an d modeling of the Galactic plane dus t con- 
tinuum: iBloemen et all 1984; 'Stro ng fc Mattoxl Il996l : 
iDame. Hartmann. fc Thaddeus 2001) . with magnetic en- 
ergy introducing a small but non- negligible correc- 
tion on the small scales (< 50%; ICrutcheJ Il999t ). 
Throughout this paper we adopt Xco = 2 x 10 ° 
cm~^(K km s~^)~^ as the Galactic conversion factor. 
Note, however, that some nearby lines-of-sight at high 
Galactic latitudes may exhibit significant larger values of 
this coefficient (e.g..lMagnani. Blitz, fc Wouterlootlll988l : 
iGrenier. Casandjian. fc Terrieill2005l )^ 

Therefore, mass estimates of GMCs show that they 
are gravitationally bound and in approximate virial equi- 
librium, which requires M^ir ^ Ra^ oc Mium- As a 
consequence of Equation [T] and virial equilibrium, all 
GMCs have approximately the same mean column den- 
sity: iVn ~ 1.5x 10^^ cm~^, or an equivalent mass surface 
density E « 170 Mq pc~^ (after applying a 40% correc- 
tion for the mass contribution of elements heavier than 
Hydrogen). An equivalent result is that all Milky Way 
GMCs have approximately the same mean visual extinc- 
tion from edge t o edge, Ay ^ 7.5, u sing the standard 
d ust-to-g a s rati o (|Bohlin. Savage, fcl ^rakc 1978). 

iMcKed ([1999) argues that Equation |4] is a consequence 
of the interplay between CO chemistry, the Milky Way 
radiation field, the dust-to-gas ratio, and the local pres- 
sure of the ISM. As a result, H2 clouds observable in CO 
are at least marginally bound. Furthermore, McKee ar- 
gues that the constancy of iVn stems from the fact that 
clouds with lower column density would not be bright 
in CO, while clouds with higher column density would 
rapidly collapse. 

Another consequence of the Larson relations is that 
Galactic GMCs have low volume-averaged densities, 
50 < n(Il2) < 500 cm~^. Therefore their dumpiness 
must be high in order to excite the rotational transi- 
tions of CO and of other commonly observed high-density 
tracer molecules such as IICO+, which require volume 
densities in excess of 10* cni^"^. This dum piness is prob- 
ably induced by supersonic turbulence (jPadoan et al.l 
119981 lOstriker. Gammie. fc Stond [19990 . Further, the 
average density of a GMC decreases with increasing 
mass, and as a r esult the free fall tim e in th e cloud in- 
creases. This led iKrumholz fc McKed ()2005[ ) to predict 
that larger GMCs will form collapsing cores and stars 
less rapidly than smaller ones. 

The purpose of this study is to determine, using a 
consistent methodology, the relations between size, ve- 
locity dispersion, and luminosity for extragalactic GMCs 
that are encapsulated in the Milky Way by the Larson 
laws. We do so by studying GMCs in dwarf galaxies and 
comparing them with those in spiral galaxies. Dwarf 
galaxies are both numerous and nearby, and host physi- 
cal conditions that fall outside those explored by surveys 
of the Milky Way. The range of parameter space thus 
probed provides an ideal laboratory to test theories of 



cloud structure and formation. Particularly, dwarf galax- 
ies have low metallicities and their lack of internal extinc- 
tion results in intense radiation fields, reminiscent of the 
conditions in primitive galaxies. They display irregular 
morphologies and are not dominated by spiral structure. 
Dwarfs often have slowly rising, nearly solid-body, ro- 
tation curves; this translates into small rotational shear 
and perhaps changes the sources and scales of turbulence 
in their ISM. The stellar potential wells of dwarf galaxies 
are weaker than those in large spiral galaxies and the am- 
bient ISM may be correspondingly less dense. These en- 
vironmental changes may influence the equilibrium con- 
ditions of GMCs, and affect their Larson relations. 

The integrated properties of dwarf galaxies already 
show clear evidence that environment affects at least 
the CO content. Dwarf irregular galaxies are fainter in 
CO than large spiral galaxies. Galaxies with metallicity 
12 -|- logO/H < 8.0 are seldom detected in CO emission 
(jTavlor et al.| [l998') . Dwarf galaxies of higher metallicity 
detected in CO, however, show similar CO luminosity 
relative to eith er their stellar or atomic gas content as 
larger galaxies (jLerov et al.ll2005l ). Here we investigate 
whether the resolved properties of GMCs in these galax- 
ies also reflect their changing environments. 

1.1. Resolving extragalactic molecular clouds 

The study of resolved GMC properties in other galax- 
ies remains technically challenging. A typical size for a 
Galactic GMC is ~ 40 pc, which at the ~ 4 Mpc distance 
of the M 81 group translates to an angular size of ~ 2". 
Achieving such a resolution requires an antenna 250 me- 
ters in diameter operating at A = 2.6 mm, the wave- 
length of the ground transition of ^^CO. This is clearly 
an impractical proposition, but aperture synthesis tech- 
niques used by millimeter-wave interferometers routinely 
attain this resolution. They do so, however, at the price 
of reduced surface brightness sensitivity with respect to 
their filled aperture equivalents, and studies such as the 
one presented here are ultimately limited by sensitivity. 
Until the advent of the Atacama Large Millimeter Array 
(ALMA) these observations will remain challenging. 

Because of these difficulties, the resolved properties of 
GMCs in galaxies other than the Milky Way remain a 
largely unexplored domain. Because of its proximity, lo- 
cated at a distance where even a modest telescope can 
attain a resolution of ^ 30 — 40 pc, the first galaxy 
where observations capable of clearly resolving individ- 
ual GMCs were po ssible was the Large Magellanic Cloud 
([Israel et al.lll982D . The first resolved study of a GMC 
beyond the immediate vicinity of the Milky Way was per- 
formed in M 31 (Vogel, Boulanger, & Bah 1987), while 
the first large scale interferometric sur vey of GMCs in 
another galaxy was done on M 33 (|Wilson fc Scovilli 
Il990t) . The first comprehensive analyses in the Magel- 
lanic Clouds were perform ed on the Small M agellanic 
Cloud (jRubio et al. 1991; R ubio. Leaueux. fc Boulanger 
ll993bD. 

iBlitz et all ()2007l ) recently presented a study of GMCs 
in the Local Group, analyzing observations from CO 
surveys complete down to a known mass limit over a 
significant fraction of each galaxy. This study con- 
cluded that the Larson relations in the Milky Way disk 
were approximately those observed throughout the Local 
Group, although some offsets (possibly methodological) 



were present. They also concluded that GMCs within a 
particular galaxy have roughly constant surface density, 
SgmCi and that once corrections for the local CO-to-Il2 
conversion factor are applied the sample as a whole has 
a scatter of only a factor of 2 in Sgmc- Finally, they 
calculated a typical conversion factor Xco ~ 4 x 10^° 
cm~^(K km s~^)^^, accurate to within 50% for most of 
the sample except the Small Magellanic Cloud, for which 
they obtain a value '^ 3 times larger. 

This paper proceeds as follows: in Sj^lwe describe the 
data sets and and galaxies studied. In ^we summarize 
the method used to consistently derive GMC properties 
across the different datasets. In ^we discuss the size- 
line width, luminosity-line width, and luminosity-size re- 
lations in our sample of extragalactic GMCs, and com- 
pare them with the Galaxy. In fj5]we discuss the implica- 
tions of these results on the equilibrium of clouds in the 
Small Magellanic Cloud, the photoionization-regulated 
star formation theory, the CO-to-H2 conversion factor, 
and the brightness temperature of extragalactic clouds. 
In each case we contrast our results in dwarf galaxies 
with those obtained in the Milky Way and the other Lo- 
cal Group spirals, M 31 and M 33. In S|6]we summarize 
this study and present our conclusions. 

2. DATA 

We use data from four tel escopes: the Ber keley- 
lUinois-Maryland Array (BIMA: rV\^lch et al.|[T99l) . Cal- 
tech's Owens Valley Radio Observatory (OVRO) 
limeter array, the Institut de Radioastronomie 
limetrique Plateau dc Burc Interferometer (PdBI), 
the Swedish-ESO SubmiUimeter Telescope (SEST) 
measure GMC properties in 11 dwarf galaxies: 
Large and Small Magellanic Clouds, IC 10, NGC 185, 
NGC 205, NGC 1569, NGC 2976, NGC 3077, NGC 4214, 
NGC 4449, and NGC 4605. We compare GMC properties 
derived in dwarfs to those found in the Milky Way (data 
obtained by the Five Colleges Radio Astronomy Obser- 
vatory, FCRAO; I Sanders et al.llT986[ ) and the other two 
Local Group spiral galaxies, M 31 and M 33. The lat- 
ter two act as key "control" data sets because they are 
obtained by the same telescopes at similar angular and 
spatial resolutions to the dwarf galaxy GMC data. Table 
[T] gives the source, telescope, line observed, resolutions, 
and reference for each data set. 

When we analyze observations of the CO J = 2 ^ 1 
transition, we assume /(2 ^ 1) = /(I -^ 0) (i.e., ther- 
malized optically thick emission). For the GMC complex 
in the N83 region of the Small Ma gellanic Cloud, this is 
within 10% of the observed ratio (jBolatto et all [20031) . 
In the Milky Way the observed ratio varies from uni ty to 
/(2 ^ 1) ~ 0.65/(1 -^ 0) re.g.. lSawada"eiraI][2f)0lh . 

Several data sets appear here for the first time. No- 
tably, GMCs in NGC 3077, NGC 4214, and NGC 4449 
where mapped at BIMA as part of a large project study- 
ing the molecular ISM in dwarf galaxies. The BIMA 
maps of NGC 2976 and NGC 4605 were obtained as part 
of the same proj ect and have appeared previously in kine - 
matic analyses (jBolatto et all 120021 : IShnon et al.ll2003^ . 
The observation and reduction procedure for these data 
follow the descriptions in those papers. 

Table [2] lists the properties of the galaxies that we 
study. Most have late type morphologies; the exceptions 
are NGC 185 and NGC 205, Local Group dwarf ellipti- 
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TABLE 1 
Data Sets Used In This Paper 



TABLE 2 
Galaxy Properties 



Source Instr. Transition Resol. Resol. Rof. 

(pc) (km s"-"-) 



Galaxy 



Morph. Dist. Mb Met. Ref. 



IC 10 


OVRO 


1-0 


19 


0.7 


5 


NGC 185 


BIMA 


1-0 


15 


4 


1 


NGC 205 


BIMA 


1-0 


28 


3 


1 


SMC N83 


BEST 


2-1,1-0 


11,17 


0.25 


9 


SMC LIRS36 


BEST 


2-1 


7 


0.05 


10 


SMC LIRS49 


BEST 


2-1 


7 


0.05 


10 


LMC N159 


BEST 


2-1 


6 


0.25 


11 


NGC 1569 


PdBI 


2-1,1-0 


22,45 


1.6 


8 


NGC 2976 


BIMA 


1-0 


94 


3 


2 


NGC 3077 


BIMA 


1-0 


41 


3 


3 


NGC 3077 


OVRO 


1-0 


41 


1.3 


6 


NGC 4214 


BIMA 


1-0 


80 


3 


3 


NGC 4214 


OVRO 


1-0 


64 


1.3 


7 


NGC 4449 


BIMA 


1-0 


117 


3 


3 


NGC 4605 


BIMA 


1-0 


109 


3 


4 



Disk Galaxies 



Milky Way 
M 31 
M 33 



FCRAO 
BIMA 
BIMA 



3-11 
27 
30 



0.65 
2 
2 



12 
14 
13 
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&) [ Walter et al. (2002): (7) Wa lter et all J2001f): ( 8) Taylor ctlil 
I999I): OrL Bolatto ct al. (2003) ; (10 ) IRubio et all lil993a'): fill 
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120031) : (liTlRosolowskv (2007) 



cals mapped by [Yound ()2001[ ) using BIMA. All of our 
sources are < 5 Mpc away; the most distant galaxy is 
NGC 4605 at 4.26 Mpc. We indicate the absolute blue 
magnitude of each galaxy corrected by Galactic extinc- 
tion to illustrate their luminosity. The galaxies span the 
range — 18 < Mb ^ — 15 magnitudes. We give metallic- 
ities for each galaxy, with the corresponding references 
for metallicity and distance. 



2.1. Dwarf Galaxies 

In this section we briefly summarize the sources and 
datasets used. For more in-depth descriptions we point 
t he reader to th e original references. 

lYound (|2001h used the BIMA interferometer to map 
CO J = 1 ^ emission in the Local Group dwarf el- 
liptical galaxies NGC 185 and NGC 205. These galax- 
ies, which are satellites of M 31, are unusual early- type 
dwarfs with CO emission, while the rest of our sample 
c onsists of late - type dwarf spiral and irregular galaxies. 

ISimon et al.l (|2003D presented a kinematic study based 
on BIMA observations of CO J = 1 ^ emission 
from NGC 2976, a dwarf spiral galaxy in the M 81 
group. These observations are considerably deeper than 
those from t he BIMA Survey of Nearby Galaxies (SONG; 
[Heifer et al.ll2003i ). The metallicity for this galaxy is un- 
certain, but the polycyclic aromatic hydro carbon abun- 
dance and dust-to-gas ratios derived by iDraine et al.l 
()2007l) suggests it is s imilar to Galactic. 

'Wa lter et all ()2OO20 mapped CO J = 1 -> emis- 
sion from NGC 3077, a member of the M 81 group cur- 
rently undergoing a dramatic interaction with M 81 and 



IC 10 

NGC 185 
NGC 205 
SMC 
LMC 

NGC 1569 
NGC 2976 
NGC 3077 
NGC 4214 
NGC 4449 
NGC 4605 



Irr/BCD 
dSph/dE3 
E5 
Bm 
Sm 
Irr 
Be 
Irr 
Irr 
Irr 
Be 



0.95 
0.63 
0.85 
0.061 
0.052 

2.2 
3.45 

3.9 
2.94 

3.9 
4.26 



-16.7 
-14.7 
-15.9 
-16.7 
-18.0 
-17.3 
-17.4 
-17.5 
-17.2 
-18.0 
-17.9 



8.2 
8.2 
8.6 
8.02 
8.43 
8.19 
8.7 
8.85 
8.23 
8.32 
8.69 



6,19 

1,17 

1,18 

7,20 

7,20 

2,16 

13 

11 

12 

15 

14 



Large Galaxies 



Milky Way 

M31 

M33 



SB 
Bb 
Bed 



0.008 
0.79 
0.84 



-21.4 
-21.1 
-18.9 



8.7 10 
8.7 8,21 
8.4 9,21 



McCain | [T995h : 
work. 



(2) 
See 
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Note. — Distances are in Mpc. Luminosities arc cor- 
rected for Galactic extinction. Mctallicities are expressed 
as 12 -f logO/ff. 



M 82, using the OVRO millimeter interferometer. This 
galaxy was also observed by BIMA during the same pe- 
riod, and those observations appear here for the first 
time. We present GMC property measurements from 
both data sets. We derive the metallicity for this galaxy 
using the in tegrated nebular fluxes o f the bright oxy- 
gen lines bv iMoustakas fc Kennicuttl (|2006f) . employing 
the calibration of McGaugh ("1991) as parametrized by 
[Kobulnicky, Kennicutt, & Pizagno (1999). The metal- 
licity is consis tent with that der i ved using the nitrogen 
c alibration bvlKewley fc Dopital (|2002). 

IWalter et all (|2001[ ) mapped CO J = 1 ^ emission 
from NGC 4214, a nearby dwarf irregular galaxy, using 
the OVRO millimeter interferometer. This galaxy was 
also observed by BIMA, and here we use those data for 
the first time. 

BIMA mapped CO J = 1 ^ emission from the 
nearby dwarf irregular NGC 4449. The data appear here 
f or the first tim e. 

IBolatto et all (|2002f ) presented CO J = 1 -^ observa- 
tions of NGC 4605, a nearby isolated dwarf galaxy, and 
used them to study the kinematics of this galaxy. The 
metallicity for this galaxy was derived using the same 
p rocedure as for NGC 3077. 

IWalter I ()2003[) mapped CO J = 1 ^ emission from 
most of the molecular clouds in the Local Group dwarf 
irregular IC 10 using the OVRO millimeter interferome- 
ter. T hese data were analyzed in detail by iLerov et al.l 
()2006t ). where maps of individual GMCs are presented. 
IC 10 is a Local Group irregular comparable in mass to 
the Small Magellanic Cloud but presently undergoing a 
burst of vigorous star formation. 



iTavlor et all (|1999[ ) used the PdBI to map both CO 
J = 2 — > 1 and J = 1 — > emission from the nearby 
dwarf irregular NGC 1569. This galaxy, which is about 
the mass of the Small Magellanic Cloud, experienced a 
starb urst period that e nded ~ 5 Myr ago. 

Bolatto et al.l ()2003l ) presented high signal-to-noise 
mapping of the N83 region in the Small Magellanic Cloud 
(SMC) using the SEST telescope. They observed both 
the CO J = 2 ^ 1 and J = 1 ^ transitions. N83/N84 
is a bright star forming region in the eastern wing of the 
S mall Magellan i c Clou d. 

iBolatto et aT\ ()2000f ) used SEST to map CO J = 2 -> 
1 emission in the N159 region of the Large Magellanic 
Cloud (LMC). This region lies just south of the bright 
star forming region 30 Doradus, and comprises clouds in 
different stages of evolution. 

2.2. Large Galaxies 

In this study we use GMCs properties measured in 
the spiral disk galaxies of the Local Group as compar- 
ison and control datasets. As discussed in the intro- 
duction, several studies ha ve established the scalings for 
Galactic molecular clouds ('Larson' 'lOSl"; 'Sol omon et al.l 
1987 : ElmcRTCcn & Fakaronc 1996; Hcvcr c t~all 120011 : 
Hever fc BrunU 120041 : IRosolowskv et al.l 1200^~Recent 



interferometric work has shown that the GMCs in the 
other disk galaxies in the Lo cal Group (M 31 and 
M 33) follow the same scalings (Iwllson fc Scovill dll990l: 
Rosolowskv et al.l 120031: iRosolowskvl 12007 : Shct h et all 
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The observations of GMCs in M 33 employed here 
were made using the BIMA millimeter interfe rometer 
and are presented by IRosolowskv et al.l ()2003l ). That 
study analyzed 20 fields in that galaxy with a resolution 
r^ Q" r^ 20 p c, targeting sources identified in the all- 
disk survey of Engargiola et al. (2003). The M 31 GMCs 
belong to a region along the nor thwest spiral arm of this 
galaxy, also observed by BIMA (|Rosolowskvll2007D . This 
region was chosen for its favorable observing geometry, 
strong dust lanes and active star formation. Both studies 
concluded that the GMCs in these disks were indistin- 
guishable in their aggregate properties from the clouds 
in the inner Milky Way. These observations constitute 
the most complete data sets for these galaxies, have sen- 
sitivity and resolution comparable to the dwarf galaxy 
data, and have been analyzed in an identical fashion. 

For the Mi lky Way we einploy t he catalog of CMC 
properties bv ISolomon et al.l (|1987f ). This catalog con- 
sists of size, line width, and luminosity measurements 
for 273 Galactic GMCs and constitutes a representative 
dat a set of Milky Way CMC properties. Like us (see 
^ , ISolomon et al.l use moments to measure cloud prop- 
erties and isosurface-based methods to identify individ- 
ual GMCs. Their study, however, does not make an ex- 
plicit correction for the effects of sensitivity (they use 
boundary isosurfaces of brightness temperature T^ = 4 
to 7 K). It also uses an outdated Galactocentric distance 
for the Sun (10 kpc). It is unclear how to best correct cor- 
rect their catalogued cloud properties to bring them into 
closer agreement with our methodology. We use here the 
original numbers, and rely on the M 31 and M 33 data as 
a consistency check. As is apparent in Figures (U |2l and 
121 these data are in good agreement with the standard 
Larson relations derived for the Milky Way sample. 



3. METHODOLOGY 

The challenge of any study that uses different instru- 
ments to study galaxies with a range of intrinsic prop- 
erties and distances is to control the systematic biases 
introduced by the differences in the datasets. We mea- 
sure clou d properties using th e algor ithm described in 
detail by IRosolowskv fc Lerovl (j2006f ) and summarized 
here. This method measures luminosities, line widths, 
and sizes of GMCs in a manner that minimizes the bi- 
ases introduced by signal-to-noise and resolution. These 
biases are the Achilles heel of comparisons of GMC mea- 
surements across different galaxi es obtained by different 
telescopes. IRosolowskv fc Lerovl proposed that moment 
measurements combined with beam deconvolution and 
extrapolation represent a robust way to compare hetero- 
geneous observations of molecular clouds. The algorithm 
is part of a package called CPROPS implemented in IDL 
and is available upon request. The remainder of this 
section summarizes the approach step-by-step. 

The robustness, reliability, an d consistency of the 
CPROP S methodology was tested bv lRosolowsky fc Lerovl 
()2006[ ) for a globally applicable set of algorithm param- 
eters. They found the algorithm performed very well for 
data with signal-to-noise ratios (S/N) > 10, but cau- 
tioned that measurements of the properties of GMCs ob- 
served with low S/N remain uncertain, particularly when 
using the default set of algorithm parameters. Unfor- 
tunately, CO is faint in dwarf galaxies, which makes it 
difficult to attain the optimal combination of S/N and 
spatial resolution. Therefore, only a few of the datasets 
discussed here meet the S/N> 10 criterion. Nevertheless, 
the data assembled here represent the best observations 
of GMCs in dwarf galaxies to date. As a consequence of 
the low S/N data, we made minor adjustments described 
below to the standard CPROPS parameters to improve the 
reliability of our results. 

3.1. Identify Signal In Each Data Set 

First, we identify regions of significant emission in each 
data set. For most of our data these are defined as regions 
that contain pixels which satisfy the condition of having 
3 consecutive velocity channels above 4cr significance (i.e. 
a high significance "core"). We extend these regions to 
include all adjacent data which has at least 2 consecutive 
channels above 1.5a significance. For a few data sets this 
did not yield the best results and we adjusted either the 
core or the edge conditions: a higher edge threshold was 
used in the high S/N N83 and N159 data sets; lower core 
thresholds were needed for the low S/N NGC 185 and 
NGC 205 data sets; a few other minor adjustments were 
also made. These changes yielded stable properties for 
clouds of low S/N, and suppressed diffuse emission blend- 
ing clouds together in the extremely high S/N case. The 
signal identification step, often referred to as "masking" , 
can have a large impact on derived cloud properties and is 
difficult to motivate physically, especially in the low S/N 
regime common to much of our data. Thus, by necessity 
this step introduces a certain amount of subjectivity in 
the analysis. 

3.2. Apportion Emission Into Clouds 

Once signal is identifi ed, data is apportioned into indi- 
vidual clouds following IRosolowskv fc Lerovl (|2006D . In 
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all cases, we use the "physical priors" described in their 
appendix, which set the parameters of the decomposi- 
tion to common values motivated by the properties of 
Galactic GMCs. This entails decomposing emission into 
clouds using effective spatial and velocity resolutions of 
15 pc and 2 km s~^ and ignoring even high-contrast 
substructure below these scales. When it is not possi- 
ble to achieve these resolutions because of limitations 
in the data, we use the actual spatial or velocity reso- 
lution of the observations. When two potential clouds 
share an isosurface, the algorithm adopts a conservative 
approach to separating them. The clouds are considered 
separate only if all of the following conditions are met: 1) 
each is large enough to measure meaningful properties, 2) 
each shows a minimum contrast between peak and edge, 
and 3) the decision to break the isosurface into multiple 
clouds has a significant effect on t he derived cloud prop- 
erties. iRosolowskv fc Lerovl ()2006D showed this approach 
to be conservative and robust in the presence of noise, 
minimizing the identification of low-signal, marginally re- 
solved (and thus often spurious) clouds. 

3.3. Use Moments to Measure Cloud Properties 

After identifying the individual clouds, we calculate 
their size, velocity dispersion, and CO luminosity. We 
determine the major and minor axis of the cloud us- 
ing principal component analysis and measure the second 
moments of emission distribution along these axes. The 
RMS size of the cloud (ar) is the geometric mean of these 
two moments. The velocity dispersion (cr„) is measured 
from the second moment along the velocity axis, and the 
luminosity via the zeroth moment (sum) over the cloud. 

At finite signal-to-noise the second moment will un- 
derestimate the true cloud size, introducing a sensitivity 
bias. Although for Galactic data it is possible to measure 
reliable properties by adopting a fixed brig htness bound- 
ary for clouds (e.g., ISolomon et al.|[T987[ ). it would be 
impossible to compare measurements of different galaxies 
without removing this bias. We avoid the sensitivity bias 
by measuring the size, velocity, and luminosity as a func- 
tion of intensity isosurface and extrapolating this to the 
case of infinite signal-to-noise. The size and line width of 
the cloud are extrapolated li nearly, the luminosity is ex- 
trapolated quadratically (see IRosolowskv fc Lerovll2006l . 
for a justification of the extrapolation orders). We dis- 
cussed in ii2.2l thc corrections that we apply to the Galac- 
tic sample to bring it into agreement with our method- 
ology. 

It is often the case for extragalactic observations that 
the size and line width of a GMC are comparable to the 
angular and velocity resolution of the telescope. We cor- 
rect for this by "deconvolving" both the telescope beam 
and the velocity channel profile. We do so by subtract- 
ing their values from the extrapolated moment measure- 
ments in quadrature. Clearly this step does not come 
for free in signal-to-noise. Size measurements of clouds 
that are only marginally resolved suffer from substantial 
uncertainty, and this is accounted for in the error esti- 
mates. 

3.4. Derive Physical Quantities 

From the moment measurements we convert to phys- 
ical units — radius, luminous mass, virial mass. The 
velocity dispersion is given directly by the moment. 



We convert the RMS size to the spherical radius of the 
cloud usin g the factor 3.4 / -^/7r, s o that radii in this paper 
follow the ISolomon et al.l (|1987t ) definition: R = l.9lar- 
This factor can be motivated by considering a constant 
density spherical cloud and comparing its RMS size to 
its radius. 

We compute luminous masses, Mium, using the for- 
mula 

Mi^m = l-l ^ 10^ Sco D^ Mq, (5) 

where 5'co = / Ico c^^ dv is the integrated ^^CO fiux 
of the molecular cloud measured in Jy km s~^, D is 
the distance to the source in Mpc, and the coefficient 
corresponds to our adopted conversion factor Xqo = 
2 X 10^° cm^^(K km s^^)^^, and includes the contri- 
bution of the cosmological He abundance to the total 
mass. Equivalently, for a CO luminosity Leo expressed 
in K km s~^ pc^, 
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We compute the virial masses under the assumption 
that each cloud is spherical and virialized with a density 
profile of the form p ex r^^. This assumption is made for 
consistency with previous work, and we cannot assess its 
validity with the dat a in hand. Thus the virial mass is 
given by the formula (jSolomon et al.lll987( ) 

M„, = 1040^2 i? Mo, (7) 

where cr„ is the cloud velocity dispersion in km s^^, i? is 
the spherical cloud radius in pc, and M^r is the cloud 
virial mass in Mq. Implicit in this equation is the as- 
sumption that the virial parameter equals one. 

Other works, in particular those focused on Milky 
Way clouds, have treated clouds as ellipsoids and used 
the corresponding form of the virial theorem (e.g. 
iBertoldi fc McKee[ll992| ). instead of the simplified Equa- 
tion [71 Because of the low resolution and S/N of extra- 
galactic observations, this approach is untenable. Un- 
der these conditions, deconvolution becomes very diffi- 
cult and the precise shape of the source is hard, if not 
impossible, to recover. We note that the more com- 
plex measurements o f clou d dimensions reduce to the 
iRosolowskv fc Lerovl (l2006f) RM S size parameter for use 
in the IBertoldi fc McKed (|l992D formulation of the the 
virial theorem. 

To relate the mass surface density, SgmCj to the co- 
efficient of the size-line width relation, cr„ = C B?'^ , we 
simply divide Equation |7| by the area of the cloud, it B?, 
to obtain 



Sgmc — 331 C , 

where the coefficient C is in units of km s^^ 
the surface density is expressed in M© pc^'^. 



pc 



(8) 
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3.5. Estimate Uncertainties 

The errors quoted in this paper are derived from boot- 
strapping the data in each cloud and re-deriving its prop- 
erties. This yields a realistic assessment of the uncer- 
tainty once signal has been apportioned into clouds. We 
also include a 25% gain uncertainty in the luminosity, 
typical of the uncertainty in the flux calibration at mil- 
limeter wavelengths. In high S/N cases, e.g. our N 83 



and N 159 data sets, bootstrapping yields very low un- 
certainties. In these cases, virtually all of the uncertainty 
in the derived cloud properties rests with the apportion- 
ment of emission into clouds. This uncertainty is not 
included in our results. 

4. RESULTS 
4.1. Derived Cloud Properties 

Table |3] presents our measurements of GMC proper- 
ties in dwarf galaxies (available in full in the electronic 
edition). In order the columns list: (1) the galaxy that 
contains the GMC; (2) a number identifying the cloud; 
(3) the transition, CO J = 2 -^ 1 or 1 — > 0, being ob- 
served; (4) the right ascension and (5) declination of the 
intensity-weighted cloud center; (6) the mean LSR veloc- 
ity of the cloud in km s^^; (7) the log of the cloud radius, 
in parsecs; (8) the log of the line of sight velocity disper- 
sion of the cloud, in km s^^; (9) the log of the luminous 
mass from Equation[5l (10) the log of the virial mass from 
Equation[7l and (11) the peak brightness temperature in 
the cloud, in Kelvin. Uncertainties from bootstrapping 
are listed with each measurement. 

Table |4] summarizes some of our results on a per galaxy 
basis. In order the columns list: (1) galaxy name; (2) 
number of measured GMC sizes in either CO transition; 
(3 — 4) average and range of GMC masses, in units of 10^ 
M©; (5 — 6) error- weighted average and range of the log- 
arithm of Xcoi in units of our adopted Galactic value; 
(7—8) error- weighted average and range of the coefficient 
of the size-line width relation, ayR~^'^; and (9) corre- 
sponding average GMC surface density from Equation 

m 

We checked for correlations between parameters in Ta- 
bleland distance or resolution, and found none obvious. 
The major bias present, which is unavoidable, is that 
there is a correlation between our minimum virial mass 
and distance that is simply due to sensitivity. In that 
sense we sample different ranges of GMC mass in dif- 
ferent galaxies. Considerations such as the nonlinearity 
of the mass-luminosity relation, for example, become im- 
portant in the Magellanic Clouds where the cloud masses 
are small. 

4.2. The Extragalactic Size-Line Width Relation 

Of the three Larson relations perhaps the most funda- 
mental one is that between cloud size and line width 
(Equation [T]), as it is mostly (but not entirely) inde- 
pendent of the details of the excitation of the tracer 
molecule. The initial studies already pointed out the con- 
nection between the size-line width relation and expec- 
tations from the Kolmogorov theory for subsonic turbu- 
lence in incompressible fluids (|Larsonlll979lll98l[ ). Later 
measurements in the inner Galaxy allowed a better de- 
termination of the relationship, showing that the ex- 
ponent is not 1/3 as it would be expected for a pure 
Kolmogorov cascade , but rather the stee per 0.5 ± 0.05 
used in Equation [T] (jSolomon et al.l Il987t ) . An index of 
0.5 is consistent with expectations for compressible su- 
personic magnetohydrodynamic turbulence (also known 
as Burger's turbulence, [y cstuto. Ostrikcr. & Stone 2003; 
iBrunt fc Mac Low! l2004f ). iHever fc Brunt! (|2004[ ) show 
that this relationship extends also within GMCs, which is 
consistent with the multi-scale turbulence interpretation 



of the observed line widths. The authors argue that be- 
cause the size-line width relation is well-defined it implies 
the invariance and universality of the turbulent scaling 
relations. For clouds in virial equilibrium M r^ RcTy, thus 
the observed size-line width relation implies M ^ R^, or 
equivalently that all clouds have similar surface densities. 

Recent studies of cloud properties in the outer Galaxy 
have shown that the size-line width relation breaks down 
for very small clouds. In these objects the velocity dis- 
persio n ay become s a constant independent of their size 
(Hever et al.l 120011 ). These authors find that outer disk 
clouds smaller than ~ 7 pc in radius have an almost con- 
stant velocity dispersion ay ~ 0.8 km s~^, while clouds 
larger than ^ 7 pc follow the usual scaling. iHever et al.l 
interpret this fact in terms of the dynamical equilibrium 
state of the small clouds. These clouds have a virial pa- 
rameter a — Myir/Miurm representing their ratio of ki- 
netic to gravitational energy, that is a > 1 for masses 
below M -- 1.3 X 10"* M©. Thus these small clouds 
have an excess of kinetic energy over their self-gravity, 
are not in virial equilibrium, and in order to be bound 
they have to be confined by an external pressure of order 
P/k ^ 10* cm~'^ K. The precise location of the break 
in the size-line width relation is then likely a function 
of the local pressure of the cloud environment. Simi- 
larly, analysis of the ^'^CO emission from GMCs in the 
inner 3 kpc of the Milky Way reveals a less steep size- 
line width relatio n than that obtain ed from ^^CO studies, 
ay - i?o.i5-o.30 (|Simonet al.ll2001f ). These authors con- 
clude that although cloud complexes are bound, many of 
the clumps identified in their study are not gravitation- 
ally bound, probably due in part to the meth odology of 
their clump identification algorithm. Finally, lOka et aTl 
(|2001h measured the GMC properties near the Galactic 
Center, finding that those clouds have line widths a fac- 
tor of 3 to 5 times larger than similar clouds in the disk 
and suggesting that this is due to the large pressure in 
this region. These studies show that there are observed 
departures from the traditional size-line width relation in 
the Milky Way, in the direction of larger line widths for 
a given size. 

Figure [T] shows our results for the ensemble of extra- 
galactic clouds. Remarkably, dwarf galaxy cloud mea- 
surements are broadly compatible with the Milky Way 
size-line width relation, represent ed by the gray a rea and 
the small dots from the sample of lSolomon et al.| |l987). 
Detailed inspection of the data reveals, however, that 
there is a tendency for GMCs in dwarf galaxies to fall 
preferentially under the Galactic relationship. This is 
most apparent in a few clouds in the LMC, IC 10, and 
most clouds in the SMC, which are significantly inconsis- 
tent with the Milky Way. Note that these departures oc- 
cur in the direction opposite from those described in the 
previous paragraph. While we may suspect that cloud 
measurements yielding sizes i? > 100 pc are affected by 
blending (such as those on the right side of the diagram), 
we have no reason to think that there are problems with 
the LMC, SMC, or IC 10, where we have some of the 
best spatial resolutions in the dataset. Furthermore, the 
control sample of M 31 and M 33 GMCs is entirely consis- 
tent with the Galactic size-line width relation and does 
not suffer from the same offsets. We discuss the par- 
ticular situation of the SMC, which exhib its the most 
convincing departures, in §5.41 Note that iRubio et al.l 
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Fig. 1. — Extragalactic size-line width relation. The diflFerent symbols and the associated error bars indicate the measurements of 
molecular clouds in different galaxies. The black dots correspond to the sample of Milky Way clouds by Solomon et al. (1987). The open 
and filled symbols indicate measurements based on CO (1 — 0) and (2 — 1) observations respectively. The dashed black line and the gray 
area indicate the standard relationship for the Milky Way and its uncertainty, ct^ = (0.72±0.07)-R'^'^^''''^. The dashed blue line illustrates 
the fit to the extragalactic GMC sample represented by Equation |9] The dash-dotted blue line indicates the fit to the dwarf galaxies only. 
The large color symbols indicate the median centroids for the GMCs in M 31 and M 33, our control sample. Dwarf galaxies are mostly 
compatible with the Galactic size-line width relation, with a tendency to deviate in the direction of smaller velocity dispersions for a given 
cloud size. In this regard the most discrepant points belong to the Small Magellanic Cloud and IC 10, some of the lowest metallicity 
galaxies in the sample. This is not, for the most part, a bias introduced by our methodology as the M 31 and M 33 points appear to follow 
the Galactic relationship. 



()1993br ) arrived at a different conclusion in tlieir study 
of tfie SMC. Because tlie GMC samples do not substan- 
tially overlap it is difficult to assess the source of the 
discrepancy. We attribute it mostly to the differences in 
the analysis methodologies, albeit it is possible that the 
samples are intrinsically different (perhaps due to differ- 
ent radiation environments). 

Although it is not entirely correct to derive a unique 
size-line width relation for the ensemble of extragalactic 
data, we have determined a best fit relation. In order to 
avoid driving the result with the small formal error bars 
in the Magellanic Clouds, and to better represent the 
scatter in the data, we increase the errors in quadrature 
in both variables until we obtain a best fit x^ value of 
unity. The resulting relation, illustrated by the dashed 
blue line in Figure [U is 



0.44lEJ;i«i?°-«°±"i"kms-i 



(9) 



where the data are fit in the log-log plane with simul- 
taneous errors in both variables, and the error bars are 
derived using a bootstrapping technique. If we only use 
the dwarf galaxy data, the coefRcient and the exponent 



are O.Sltj;];? and 0.65 ± 0.07 respectively. 

If the size-line width relation is interpreted in terms 
of virialized self-gravitating clouds supported mainly by 
turbulence, the compatibility between the Galactic rela- 
tion and our measurements implies that most GMCs in 
dwarf galaxies have surface densities not very different 
from Galactic GMCs, for which Sgmc « 170 M© pc~^. 
Extragalactic GMCs that fall under the Galactic size-line 
width relation would have lower surface and volume den- 
sities than corresponding clouds in the Milky Way. The 
small-size, low velocity dispersion end of the dash-dotted 
line in Figure [T] corresponds to a drop of approximately 
a factor of 2 in the velocity dispersion, which translates 
to a factor of 4 decrease in surface density with respect 
to the Galactic value. This surface density may be ap- 
propriate for the most extreme galaxies in the sample 
(c.f.. Tabled]). If we fit a model that is proportional to 
i?*^'^, the sample of dwarf galaxies is consistent with an 
average surface density Sgmc "^ 85 M0 pc~^. 

An alternative possibility is that the virial parame- 
ter, a measure of the ratio of kinetic to grav itational 
energy of the cloud (jBertoldi fc McKed Il992f ). is dif- 



ferent for these objects. The virial parameter, a, can 
be expressed in ter ms of the velocity dispers ion follow- 
ing Equation 27 in iMcKee fc Ostrikeil (|2007t ) to obtain 
a^ — VaSGMC R~°'^- Thus, a low coefficient in the size- 
line width relation can be attributed to a lower than 
Galactic virial parameter. Bound clouds in the Milky 
Way have a ~ 1, while pressure-confined clouds exhibit 
Of ^ 1 — a cloud in equilibrium with a < 1 requires 
some additional type of support, presumably provided 
by magnetic fields. On the spatial scales studied here, 
such support is usually unimportant in the Milky Way. 

4.3. The Luminosity-Line Width and Luminosity-Size 

relations 

Figures [2] and [3] show the relations between ^^CO lu- 
minosities and either velocity dispersion or cloud radius. 
These luminosities can be related to the corresponding 
luminous mass Mium for a Galactic Xqq using Eqs. [5] 
and[Sl The steep dependency on Uy in Equation[2]is a di- 
rect consequence of the observed size-line width relation 
in combination with the assumption of virial equilibrium. 
The coefficient in front of this equation reflects a combi- 
nation of the area ffiling fraction by CO clumps within 
the resolved cloud, and the kinetic temperature of the 
gas itself. 

Inspection of Figures [5] and [3] reveals that, despite the 
variety of environments sampled and the consistently 
lower-than-Galactic metallicity of our sources, there is 
good agreement between the properties of extragalactic 
and Galactic GMCs. Nonetheless, it is clear that there 
is an overall tendency for clouds in dwarf galaxies to be 
underluminous for their size, and clouds in the SMC tend 
to be also overluminous for their velocity dispersion. In 
other words, they tend to be larger and (in the case of 
the SMC) more quiescent (as traced by their velocity 
dispersion) than Milky Way GMCs of similar luminosity. 
The fit to all the extragalactic data, conducted as in 
shows that the luminosity-line width relation is 



Leo 






.3.35±0.19 K km S-1 pc2, (10) 

where the shallower-than-Galactic slope is driven mostly 
by the small SMC clouds. A similar analysis but includ- 
ing only the dwarf galaxy data yields a coefficient 688^5^7^" 
and an exponent 3.43 ± 0.22, thus very consistent with 
the results when M 31 and M 33 are included. The cor- 
responding extragalactic luminosity- size relation is 



Leo « 7.8t^:?i?2-^4±°-20 K km s'^ pc^. 



(11) 



where the corresponding coefficient and exponent when 
fitting only the dwarf galaxy data are G.Ola'g and 2.47 ± 
0.17. This shows that the data are completely consistent 
with the Galactic slope, but the relation is noticeably 
displaced toward lower luminosities for a given size. 

Figure [3] suggests that, despite our efforts to correct 
for them, there may be a remaining methodological off- 
se t in the cloud size de termination between the sample 
of ISolomon et al.l ()1987[ ) and this study, of order 30% to- 
ward larger R (equivalent to a factor of 0.5 toward lower 
luminosities). Unfortunately, establishing the origin of 
these offset would require a full reanalysis of the original 
Galactic survey. With the quality of the data at hand it 



is difficult to assert whether methodological or real off- 
sets are driving the observed shift toward larger cloud 
sizes, although it seems that the likely methodological 
offsets are too small to explain the level of discrepancy 
in many sources — most notably the SMC. 

Figure 2] illustrates the behavior of luminosity and 
virial mass described by equation H] for our sample of 
extragalactic GMCs. ISolomon et al. I discuss this empir- 
ical relation in some detail, showing that it is a nat- 
ural outcome of the combination of the size-line width 
relation with the assumption of virial equilibrium when 
all clouds share similar kinetic temperatures and area 
filling fractions in their resolved CO emission. Despite 
the aforementioned tendency for clouds in dwarf galax- 
ies to be underluminous for their sizes, their virial masses 
are, on average, only slightly underpredicted by their 
CO luminosity. In the particular case of GMCs in the 
SMC, because the observed departures from the stan- 
dard luminosity-line width and luminosity-size relations 
occur in opposite directions, their larger sizes and smaller 
velocity dispersions for a given luminosity approximately 
cancel out in the virial mass calculation. Thus most re- 
liable extragalactic clouds in our sample are remarkably 
compatible with equation 21 The fit to the entire extra- 
galactic sample yields 

while the fit to the dwarf galaxies only has a coefficient 
of 10.3t^;?, and an exponent 1.00 ± 0.05. Thus the ex- 
tragalactic and Galactic mass-luminosity relations are in 
excellent agreement for GMCs in the mass range 10"^— 10^ 
Mq. This suggests that, by comparison with Equation[Sl 
Xco is on average a factor of 1.7 to 2.3 larger than our 
adopted Galactic value of 2 x lO^*^ cm^^(K km s^^)^^ 
(although with large error bars). We will discuss the 
consequences for Xco and the presence of any metallic- 
ity trends in the CO-to-H2 conversion factor in ^ ^5.21 

4.4. Comparison with Complete Samples 

We outlined in ^ 31. II some of th e conclusions of th e re- 
cent study of the Local Group bv lBlitz et"aI1 ()2007D . We 
employed here the same datasets for M 31 and M 33, an- 
alyzing them with a lat er ver sion of the same CPROPS 
algorithm. While Blitz et al.l emphasized completion to 
study the cloud statistics, the present analysis made use 
of the available data with the best possible combination 
of sensitivity and resolution, thus employing different 
datasets for the Magellanic Clouds and IC 10. Further- 
more, this study incorporates a large number of galaxies 
in and beyond the L ocal Group that were not analyzed 
in lBlitz et al.l (|2007f ). Many of the results, however, are 
reassuringly similar: extragalactic GMCs broadly share 
the same properties observed in Milky Way GMCs, al- 
though the clouds in some galaxies (notably the SMC) 
appear to be somewhat different from the remainder of 
the sample. We will argue later ( !j5.4p that this is prob- 
ably caused by the low metallicity of the parent galaxy. 

5. DISCUSSION 

5.1. Implications for Photoionization- Regulated Star 
Formation 

iMcKed ()1989D discusses a self-regulated theory of star 
formation, where photoionization caused by the inter- 
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Fig. 2. — Extragalactic luminosity-line width relation. As in Fig. [T] the different symbols and the associated error bars indicate the 
measurements of molecular clouds in different galaxies, with open and filled symbols indicating measurements based on CO (1 — 0) and 
(2 — 1) observations respectively. The black dots correspond to the sample of Milky Way GMCs by Solomon et al. (1983)i s-^d the gray 
regions to their la dispersion around the dashed lines, which follow Equations [2] and |3] The dashed blue line illustrate the fits to the 
extragalactic GMCs. The dash-dotted blue line corresponds to the results of the fits for the dwarf galaxy data only. The large color symbols 
indicate the median ccntroids of the clouds that have R measurements for our control sample of GMCs in M 31 and M 33. 



TABLE 3 
Radio Observations 



Galaxy 


# 


Trans. 


R.A. 


Declination 


Vlsr 


logio R 


logio o" 


logio A^ium 


logio ^vir 


Tb 










{J2000) 


(J2000) 


(kms-i) 


(pc) 


(km s-i) 


(Mo) 


(Mq) 


(K) 




BIMA 


NGC185 


1 


1 — 





00''38""56.4= 


48°20'19.9" 


-194.1 




0.69 ±0.16 


4.52 ±0.18 




1.2 


NGC205 


1 


1 — 





00''40™24.0= 


41°41'52.1" 


-248.9 


1.43 ± 0.09 


0.78 ±0.11 


5.02 ±0.12 


6.01 ±0.24 


1.0 


NGC2976 


1 


1 — 





09''47™23.6= 


67°54'40.1" 


-31.3 


1.77 ±0.17 


0.84 ±0.14 


5.28 ±0.21 


6.46 ±0.30 


0.3 


NGC2976 


2 


1 — 





09'^47™24.1= 


67°54'33.9" 


-24.8 




0.25 ±0.33 


4.65 ±0.33 




0.3 


NGC2976 


3 







091147^1498 


67°54'44.0" 


-19.5 


1.75 ±0.18 


0.52 ±0.14 


5.43 ±0.11 


5.82 ±0.30 


0.4 


NGC2976 


4 







09h47""18.4= 


67°54'48.5" 


-6.2 


2.06 ±0.22 


0.82 ±0.17 


5.54 ±0.16 


6.71 ±0.28 


0.3 


NGC2976 


5 







09!^ 47"" 17. 3= 


67°54'59.2" 


5.5 


1.96 ±0.25 


0.29 ±0.30 


5.49 ±0.30 


5.56 ±0.51 


0.3 


NGC2976 


6 







09'"47""15.3= 


67°55'04.9" 


15.4 


2.09 ±0.19 


0.63 ±0.18 


5.88 ±0.23 


6.37 ±0.28 


0.4 


NGC2976 


7 


1 — 





09ii47mi53s 


67°55'16.7" 


20.8 




0.68 ±0.29 


5.51 ±0.36 




0.3 


NGC2976 


8 


1 — 





09'"47'"13.0= 


67°54'54.5" 


-2.6 




0.63 ±0.13 


5.18 ±0.20 




0.3 



Note. — Table [3] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown he 
guidance regarding its form and content. 
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Fig. 3. — Extragalactic luminosity-size relation. Symbols and labels as in Fig. |2] 



TABLE 4 
Galaxy-wide Measured GMC Properties 



Galaxy 


Meas. 


M 


.>r/10'^ 




log [Xco/2 


X 1020" 
s-)-i) 


CTvR 


-0.5 




Sgmc 








(M©) 




(cm 2(K km 


(km s-i 


pc-0-^ 


) 


(Mq pc-2) 






mean 


range 


mean^ 


range^ 


moan^ 


range 




IC 10 


9 


2.6 


0.2- 


5.6 


-f0.28±0.07 


0.5-5.9 


0.57 + 0.08 


0.34- 


-0.94 


108+31 


NGC 205 


1 


10.4 






-f 0.99 ± 0.27 




1.16 + 0.25 






447+191 


SMC 


11 


0.1 


0.02- 


-1.2 


+0.43 ± 0.03 


1.1-6.6 


0.37 + 0.08 


0.18- 


-0.75 


45+21 


LMC 


7 


1.3 


0.1- 


3.1 


-f0.30±0.04 


1.1-3.0 


0.61 + 0.10 


0.35- 


-0.70 


123+39 


NGC 1569 


1 


2.1 






-1-0.52 ± 0.78 




0.72 + 0.44 






170+209 


NGC 2976 


7 


26.4 


3.6- 


-53 


+0.67 ±0.11 


1.2- 15.1 


0.57 + 0.06 


0.20- 


-0.90 


106+21 


NGC 3077 


4 


16.0 


1.9- 


-48 


+0.06 + 0.12 


0.7-1.5 


0.71 + 0.08 


0.47- 


-0.99 


166+39 


NGC 4214 


3 


6.0 


4.3- 


6.9 


-0.07 + 0.17 


0.7-1.6 


0.56 + 0.12 


0.48- 


-0.94 


102+44 


NGC 4449 


1 


85.3 






+0.08 + 0.19 




0.56 + 0.09 






105+33 


NGC 4605 


3 


24.2 


19- 


30 


+0.23 + 0.21 


1.2-2.2 


0.37 + 0.08 


0.29- 


-0.48 


44+18 


M 31 


44 


5.3 


0.6- 


-21 


+0.26 + 0.03 


0.5-5.0 


0.72 + 0.03 


0.35- 


-1.62 


172+15 


M33 


14 


2.3 


0.9- 


8.0 


+0.12 + 0.06 


0.5-2.6 


0.70 + 0.07 


0.46- 


-2.40 


160+31 



Error-weighted mean and corresponding error. Range for Xco/2 X 10" . 
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Fig. 4. — The relationship between CO luminosity, Lqo, and virial mass, M^i^, for Galactic and extragalactic GMCs. As in previous 
figures white and black symbols indicate meas urements based on CO 1 — and 2 — 1 respectively. The dots correspond to the sample 
of Milky Way GMCs by ISolomon et all I I1987IV The dashed black line illustrates equation |4] with the gray region corresponding to the 
Galactic la dispersion around the line (~ 50%). The dashed and dash-dotted blue lines illustrate the fits for all extragalactic GMCs and 
for the dwarf galaxy sample alone, respectively. 



stellar radiation field controls the rate at which clouds 
contract and form new stars. In this theory clouds un- 
dergo contraction until the energy input by the newly 
formed stars balances the gravitational collapse. Sup- 
port against contraction at the level of individual clumps 
is provided by their magnetic field. This contraction hap- 
pens at the ambipolar diffusion rate, as the neutrals dif- 
fuse through the ions which are anchored by the mag- 
netic field. The rate of diffusion is determined by the 
ionization fraction of the cloud, which in its outer regions 
(which comprise most of its mass) is set by the ultravi- 
olet interstellar radiation field (ISRF), and in its inner 
regions is due to cosmic ray ionization. The total rate of 
star formation of a cloud thus depends critically on its 
extinction and radiation environment. Clumps immersed 
in weaker ISRFs require less extinction to collapse, and 
vice versa. 

These ideas have considerable implications for the 
structure of GMCs as a function of environment, and 
it is of great inter est to t e st the ir predictions. In the 
theory outlined by iMcKed ()1989[ ) , extinction is the key 
parameter that regulates the equili brium o f clouds, not 
total column density. Accordingly, iMcKed proposed as 
one of the key tests of the theory that clouds in galaxies 
with different metallicities and dust-to-gas ratios should 
obey different size-line width relations, as 



0.72 



A, 



7Ma 



0.5 



i?"-^ km s- 



(13) 



This expression is the analog of equation [Tl where the 
numerical coefRcient now incorporates the cloud mean 
extinction. Ay, and the ratio of Ay extinction per hy- 
drogen nucleus relative to the Milky Way, 5gr , so that for 
Galactic clouds we recover the standard size-line width 
relation. The prediction from photoionization-regulated 
star formation theory is that equilibrium molecular 
clouds in galaxies with low values of Sgr should have 
high velocity dispersions for a given size. In other words, 
clouds in low metallicity environments need to be more 
massive to attain the characteristic extinction in their 
inner regions for a given ISRF. 

Figure [5] shows the average coefficient of the size-line 
width relation for each galaxy in our sample plotted as a 
function of metallicity, which is a good proxy for Sgr in 
equation [131 Neither the placement of individual GMCs 
in the size-line width relation (Fig. [T]), nor Fig. [5] show 
evidence of the systematic behavior predicted by the the- 
ory. If anything, clouds in IC 10 {6gr ~ 0.3) and the SMC 
{Sgr < 0.2), which should have higher-than-Galactic ve- 
locity dispersions by factors of ^ 2 — 3, appear to have 
systematically lower-than-Galactic values of velocity dis- 
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persion for a given size. Even if we discard the clouds 
that fall under the Galactic size-line width relation as 
potentially out of equilibrium or perhaps magnetically 
supported, the remaining clouds do not show larger than 
Galactic velocity dispersion. It is important to real- 
ize that several of these GMCs e xist in environrn ents 
of larger-than-Galactic ISRF (e.g. Jlsrael et al.|[l996[ ). so 
that the observed lack of a trend in (Tt,i?~°'^ is not di- 
rectly caused by a cancellation between the effects of 
ISRF and metallicity. 

Early support for photoionization r egulated star for- 
mation was found bv lPak et al.l (|1998( l. who used IRAS 
60 and 100 /xm data, a model of PDR regions, and H2 and 
[CII] measurements to derive the structure of molecular 
clouds in the Magellanic Clouds. Their determination of 
column densities, however, hinges critically on the esti- 
mate of the radiation field, which is derived from IRAS 
measurements. In the SMC they model dust tempera- 
tures in the r ange of 42 — 49 K, higher than the temper- 
atures in the iStanimirovic et all (2000) map which em- 
ploys the same data. We use the IRIS data to check the 
temperature towards LIRS36, one of their target regions, 
and find i t to b e 33 K instead of the 45 K derived by 
iPak et all (|1998f) . Because the estimated radiation field 
is very strongly dependent on the assumed dust temper- 
ature, this alone suggests a '--^ 5 times lower field than 
they adopt. This is enough to bring their result into 
approximate agreement with our own, and makes the ra- 
diation field in this region intermediate between Galactic 
and 30 Doradus in the LMC. A further concern with this 
analysis stems from the use of 60 /xm data to establish 
a temperature, as it is known that there is an important 
contribution from stochastically heated small grains at 
this wavelength (pa. rticularlv in the Magellanic Clouds, 
[Bernard et al.l l2008l). Studies c ombinin g 100 iim data 



with 160 - 170 /xm (JLerov et al.l [2007: Wi lke et al.l[2004D 
determine typical temperatures of 20 — 23 K for these 
SMC regions. 

While our analysis does not directly support 
photoionization-regulated star formation theory, there 
are several caveats that need to be considered before 
rejecting it. First, strictly speaking the increased col- 
umn density prediction from photoionization-regulated 
star formation theory applies to the structures that form 
stars. Thus, it is possible that these observations do not 
sample the scales at which the column density enhance- 
ments occur. Second, in addition to assuming clouds in 
virial equilibrium, this prediction assumes that the ratio 
of mean clump to cloud extinction is similar to Galactic 
clouds, and that the characteristic density (at which cos- 
mic ray ionization is « 10~^) is a similar fraction of the 
cloud density as in Galactic clouds. Unfortunately, these 
parameters are almost completely unconstrained for ex- 
tragalactic GMCs. Third, metallicity may not be a good 
linear proxy for Sgr- This parameter folds in the effects 
of dust-to-gas ratio and extinction, and although dust- 
to-gas ratio likely tracks metallicity the dependence of 
extinction on metallicity may be important. In the SMC 
N(R l)/Av » 13.2 X 10^1 cm-2 mag-\ a factor of 6.5 
larger than iMcKed 's adopted Galactic value and com- 
patible with the rati o of metallicities be tween the SMC 
and the Milky Way (I Gordon et a l."2003"l. At ultraviolet 
wavelengths, however. I Gordonet al. find that Ax/ Ay is 
larger than Galactic by factors of 2—3. As a result of this. 



and assuming that the ultraviolet extinction and not Ay 
is the parameter that matters at setting the cloud ion- 
ization level, we may probe a smaller range of Sgr than 
that assumed in Fig. [5l Note that although this consid- 
eration diminishes the discrepancy with the theoretical 
predictions it does not eliminate it. Fourth, equation 
[13] requires that CO is present throughout the molec- 
ular material so that its velocity dispersion accurately 
samples the potential of the GMC. This may not be the 
case in low metallicity environments if there are large II2 
envelopes faint in CO surrounding the CO-bright por- 
tions of the GMC, and it may constitute the biggest 
limitation in our test of photoionization-regulated star 
formation. Nonetheless, the result stands that there is 
no observed trend of increasing velocity dispersion in the 
size-line width relation for decreasing metallicity on the 
scales sampled — clouds show strikingly similar proper- 
ties in all galaxies. 

5.2. The CO-to-H2 Conversion Factor 
5.2.1. Background 

Using ^^CO observations to measure the amount of 
molecular hydrogen requires assuming a CO-to-H2 con- 
version factor Xqo = iV(Il2)/ / /(CO) dv (jLebrun et al.l 
1983) . This factor incorporates the effects of abundance, 
excitation, and cloud structure averaged over a large 
area. The optically thick ^^CO 1 — transition can 
be used to trace the molecular mass of a cloud because 
it arises from the surfaces of clump ensembles within 
the telescope beam: by measuring Ico we are counting 
clumps, and under very general conditions the CO in- 
tensity will be proportional to the total molecular mass 
of the clump ensemble. T his is the essence of th e CO 
"mist" model discussed by iDickman et al.l ()1986[) . Be- 
cause encompassing all the relevant physics of the prob- 
lem into ab initio calculations is extremely difficult and 
not well constrained, astronomers rely on empirical cal- 
ibrat ions of ^co obtained using a variety of techniques 
(e.g., Dickman 1978; Sanders et al. 1984: Bloe men et al.l 
11986; Strong & Mattox 1996; Dame et al. 2001^^ 

It is reasonable to expect that Xqo will depend on the 
local properties of the interstellar medium (ISM); vol- 
ume density, temperature, radiation field, and metallicity 
(e^ lKutner fc Leung|[l985t lElmegreenlll989l : iBell et al.l 
|2006[) . Observations throughout the Galactic disk, how- 
ever, constrain Xqo within a narrow range Xco ~ 
(1.8 ± 0.3) X 10^° cm-2(K km s-i)"\ with excursions 
of up to a fac tor of ^ 2 oyer th is value, particularly at 
high latitudes (jDame et al.ll200l[) . This suggests that the 
local variations due to volume density, temperature, and 
radiation field are unimportant when cloud properties are 
averaged over a large area, which is the case for most ex- 
tragalactic observations. In cases where the average con- 
ditions are very different from those in the Galaxy, how- 
ever, differences in Xqq become apparent. For example, 
mass estimates of molecular clouds subjected to extreme 
radiation fields, such as those in starburst environments, 
suggest that they have an Xco convers ion factor 4 to 2 
times smaller than in the Galaxy fe.g.. lYao et al.ll2003f ). 
This can be understood in terms of a larger emissivity 
per CO molecule in starburst galaxies, stemming from a 
higher excitation state due to larger physical tempera- 
tures and volume densities in these environments (e.g.. 
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Fig. 5. — CoefBcient of the size-line width relation for the different galaxies in our sample, as a function of metallicity. Each symbol 
represents the error-weighted average value of a^,/ B?-^ for all clouds in a galaxy. The horizontal bars indicate the uncertainty in galaxy 
metallicity, while the vertical fill bar correspond to the error in the error- weighted average and the dotted vertical line illustrate the full range 
of valu es for av/BP'^ for individual clouds within a galaxy. The thick lines indicate the behavior for the Galactic sample of Solomon et al] 
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the normalized extinction per hydrogen nucleus (Sgr) is proportional to the source metallicity, while the dotted horizontal line and gray 
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E21 



IWeifi etaLl 120011 ). Another possibility is that a key as- 
sumption underlying the CO-to-H2 proportionality, that 
the CO line width trac es the H2 mass, bre aks down in 
very CO-rich systems (jDownes et al.|[l993l) . Similarly, 
Xqo conversion factors a few times lower than that of 
the Galactic disk are commonly obs erved in the centers of 
spiral galaxies, including our own (jSodroski et al.|[l995t 
llsrael. Tilanus. fc Baas"2006[). 

A property of the ISM that is expected to have a dra- 
matic impa ct on the value of Xgo is its heavy-elemen t 
abundance (jlsrael et all 119861: iMalonev fc Blackl I1988D . 
Unlike density or radiation field fluctuations, metallic- 
ity does not average out over the large areas sampled 
by a telescope beam. In fact, dwarf galaxies, which 
tend to have subsolar metallicities, show little evidence 
for metallicity gradients at least in the ex amples where 
such gradients could be well ineasured (e.g..lDufourlll984l : 
iKobulnickv fc SkillmanllT997t ISkillmanl l2003D. Metalhc- 
ity affects cloud structure in two ways: 1) directly, as 
smaller abundances of C and O translate into lower CO 
yields, and 2) indirectly, as the dust-to-gas ratio is low- 
ered, which in turn diminishes the H2 formation rates 
and the shielding of molecular gas from the photodisso- 



ciating effects of ultraviolet radiation. 

The observational effects expected from these changes 
can be summarized as follows: as metallicity decreases, 
decreasing dust shielding pushes the C+ to CO transition 
further into the molecular gas. In terms of a spherical 
clump, the radius of the t = 1 surface for CO becomes 
smaller as the metallicity decreases. Because at high- 
enough extinction most of the carbon in the gas phase 
turns into ^■^CO, which in turn becomes rapidly optically 
thick, the actual C/H ratio of the gas has only a small 
effect on the position of r = 1 surface. For a typical 
Galactic cloud, for example, the C+/CO transition oc- 
curs at Av ~ 2 while the t = 1 surface of CO occurs 
at column densities A^co ~ 8 x 10^^ cm~^ (km s~^)~^ 
which translate into visual extinction of Ay ~ 0.015 for 
a typical gas-phase carbon abundance. Thus, almost im- 
mediately after becoming the dominant form of carbon, 
^^CO turns optically thick. Because of its large abun- 
dance, H2 is self-shielding in Milky Way mole cular clouds 
(lAbs-rall et al.l[l99llD7sdiie fc Bertoldil[l99fih . CO, how- 
ever, is only mildly self-shielding, although it also cross- 
shield s with H2 at column dens ities A^(H2) > 10^^ 
cm~^ (jvan Dishoeck fc Blackl[l98a l . If the conditions are 
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such that H2 self-shielding dominates over dust shielding 
at lower metallicities, the atomic to molecular H I/H2 
transition will not move into the cloud for decreasing 
dust-to-gas ratios. As a result, a molecular clump will 
not change its size with decreasing metallicity but its 
^^CO emitting core will diminish, yielding an in creas- 
ing Xqo conversion factor (e.g., Maloncy & Blacki JToSSt 
iBolatto. Jackson, fc Ingallslll999i : iR5Uig et al.ii2006l l 

5.2.2. Existing Calibrations 

There are several observational calibrations of Xco 
with metallicity, Z, in the literature, showing a range 
of behaviors. Most of the calibrations find an increas- 
ing Xco with decreasing Z although the rate of increase 
varies greatly depending on the technique and the spatial 
resolution of the observations used. 

A key datum in these calibrations is the measurement 
in the Small Magellanic Cloud, as this is the galaxy 
with the lowest metallicity where CO is reliably detected 
[ZsMC ~ 0.2Zg: lDufoinlll984t IVermeii fc van der Hulstl 
2002) . In t he firs t systematic study of this galaxy. 



Rubio et al.1 ()1993bD used single-dish CO observations 



and the assumption of clouds in virial equilibrium to 
determine that Xco ~ 60 x 10^° cm-2(K km s^i)"! 
at a resolution of ~ 160 pc, while finding Xqq ^ 
9 X 10^° cni ~ ^(K k m s^^)^^ on scales of ~ 15 pc. 
iMizuno et al.1 ()2001f ) revisited the estimate of Xqo in 
the SMC with new single-dish observations at a res- 
olution of ~ 50 pc and found Xqo ^ 10 — 50 x 



10^° cm-2(K km s' 



in rough agreement with the 



previous results. One important issue with these deter- 
minations is that it is unclear whether the large scale 
structures observed in CO are self-gravitating, virialized 
molecular clouds. The applicability of the virial theorem 
on scales of 100 pc and larger and the ability to accu- 
rately measure the size of the CO structures, both neces- 
sary to apply equation[7l are problematic. Another issue 
with these determinations is that they do not account 
for the finite angular resolution of the telescope, which 
is in many cases comparable to the sizes of the clouds. 
This will systematically bias these measurements in the 
direction of lar ger virial masses an d larger Xco- 

Subsequentlv lRubio et al.l (|2004D used millimeter-wave 
continuum observations to obtain the mass of the quies- 
cent cloud SMCBl-1 and compare it with its CO emis- 
sion, reaching the conclusion that the dust-derived mass 
is 7 to 20 times larger than the cloud virial mass. They 
suggest that this H2 is located in an envelope of the 
cloud that does no t emit brightly in CO. Very recently, 
iLerov et al.l ()2007[ ) used new far-infrared (FIR) images 
of the SMC obtained by the Spitzer Space Telescope to 
model the dust emission and locally calibrate its dust-to- 
gas ratio, obtaining a map of II2 in this galaxy. Compari- 
son with the CO emission indicates that the II2 clouds are 
on average ^ 30% more extended than the CO-emitting 
regions, so that CO clouds are immersed in an extended 
molecular envelope of II2. Furthermore, over the vol- 
ume occupied by C CLerov et al.. find Xco ^ 60 x 10^" 
cm~^(K km s^^)^^, while the overall Xqo for the entire 
SMC is approximately twice that. 

Among the calibrations for Xqq with metallicity, a 
number of them rely on the assumption of virial equilib- 
rium to obtain the molecular mass of GMCs. Wilson 
(|1995l) used interferometric CO observations of M 33 



and d warf galaxies i n conju nction with the aforemen- 
tioned I Rubio et al.l (|1993bf ) results to establish that 
Xco ~ Z~°-^'^, where Z is measured using the oxy- 
gen abundance, as O/H. lArimoto et al.1 ()1996[ ) used sev- 
eral observations in the literature to determine that 
Xco ~ Z~^. Employing some of the same interfero - 
metric observations analyzed here, IWalter et al.1 ()200lD 
determined a Galactic value for Xqo in NGC 4214, a 
galaxy where Z ^ 0.3 Zq. A result of 60% of the Galac - 
tic value was obtained in NGC 3077 ([Walter et al.ll2002f) . 
a gala xy with approximately G alactic metallicity. Sim- 
ilarly, iRosolowskv et al.l (|2003f l used interferometric ob- 
servations of the entire disk of M 33 and found no de- 
pendence of Xco on metallicity over a range of 0.8 dex 
in Z (a factor of 6). All these studies are based on virial 
mass techniques. 

Calibrations tend to be even more discrepant when 
other methods are employed. Be s ides th e afor emen- 
tioned results bv iLerov et all ()2007| ). llsraen ()1997f ) used 
H I and Infrared Astronomy Satellite (IRAS) observa- 
tions of several galaxies, estimating A"co ~ 120 x 
10^° cm-2(K km s"i)"Mn the SMC and finding a metal- 
Hcity dependence Xqq ~ Z~^'^, or Xco '^ Z~^-^ 
when taking into account the lo cal interstellar radi- 
ation field. iMadden et al.l (|1997f ) used measurements 
of the FIR [CII] (^P3/2 ^^Pi/2) transition to esti- 
mate Xco ^ 50 X 10^° in some regions of the low 
metallicity dwarf galaxy IC 10, and over 100 times the 
Galac tic value overall {Zicw ^ 0.3 Zq). Ilmara fc Blit J 
(J2007t ) used stellar extinction in the LMC to mea- 
sure Xco '^ 9.3 X 10^° cm-2(j<^ Ij-j^ s~^)~\ toward the 
low e nd of virial mass estimates for th is source. Fi- 
nally, iBoselli. Lequeux. fc Gavazzil 1)2002') found a milder 
Xco ~ Z~^ dependence using a combination of virial 
and dust-continuum methods. In summary, there are 
large discrepancies between the different authors and 
techniques, and in general estimates based on FIR ob- 
servations find stronger dependencies on metallicity than 
those based on virial arguments, although typically they 
also probe larger scales. 

A problem with the available studies has been the lack 
of uniformity in the datasets, methods, and analysis tech- 
niques. Unexpectedl y, the most uniform cafi bration in 
the literature (that of IRosolowskv et al.ll2003l ). where all 
the clouds are at the same distance, the spatial resolution 
is good (20 pc), and identical analysis is applied to the 
data finds that Xco is independent of Z. Very recent ob- 
servations, however, cast doubts o n the magnitude of the 
metallicity gradient in this galaxy. IRosolowskv fc SimonI 
()2007| ) find that the metallicity gradient of M 33 is a fac- 
tor of 3 shallower than previously accepted, which implies 
that the Rosolowsky et al. (2003) CMC data probes a 
considerably smaller range of Z than previously thought. 
Thus mild metallicity dependencies of Xco within this 
galaxy may be masked by the internal scatter of the mea- 
surements. Another potential issue is that similar gradi- 
ents in the radiation field and metallicity with galacto- 
centric distance have opposite effects and may conspi re 
to produce an almost constant Xco l|Elmegreenlll989t) . 

5.2.3. Our Results 

The fact that the overwhelming majority of our ex- 
tragalactic GMCs with reliable property determinations 
are compatible with the empirical Galactic relationship 
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Fig. 6. — Ratio of virial to luminous mass (i.e., Xco factor including the contribution from He) as a function of metallicity for the 
different galaxies in the sample. The ratio is normalized to the Galactic Xqq, as all luminous masses are computed using our adopted Milky 
Way value of the CO-to-H2 conversion factor {Xqq = 2 x 10^" cm~^(K km s~^)~^, shown by the horizontal dotted line). The horizontal 
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averages in each galaxy, and the dotted lines show the range of ratios found for all clouds in each galaxy with more than one identified 
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The thick lines show the results for the ISolomon et al.l l|1987l) sample of Galactic clouds with masses larger than 10^ Mq. The dashed 
line illustrates Xcq{Z) cx Z~^ . Note that the departure from the standard Galactic Xqq for the SMC is in accordance with the Galactic 
relation between luminosity and virial mass (see Equation|3J. Thus it is not directly associated with metallicity, just a consequence of the 
smaller average mass of the SMC clouds: similar mass clouds in the inner Milky Way show the same ratio of virial to luminous mass. 



between Lqq and M^^^ described by equation 2] fi j4.3p 
implies that we do not observe extreme departures from 
a Galactic Xqo- Indeed, Figure [6] shows our results for 
the extragalactic calibration of Xco with metallicity, us- 
ing virial GMC masses. The gray region illustrates the 
approximate range of Xco in the Milky Way found by 

1 Dame et al.l (|2001[ ). a factor of 2 around the nominal 

2 X 10^° cm-2(K km s-^)-^ The thick lines show the 
disper sion of values for the GMCs in the ISolomon et al.l 
(|1987f ) sample with masses larger than 10^ M0. The 
symbols show the average of Mmr/Mium for all GMCs in 
each galaxy in our sample, with the corresponding ver- 
tical bars indicating the error and range of Myir / Mium 
within a galaxy. The galaxy that shows the largest de- 
parture from Galactic Xco is NGC 205, but with only 
one identified cloud we lack significant statistics. The 
situation is similar for NGC 1569. Note that although 
the GMCs in the SMC show an average Xco three times 
Galactic, they are not unusual compared with Galactic 
clouds of the same mass (recall that Equation |4] is not a 
linear relation, and Galactic clouds with masses similar 
to our SMC clouds have similar ratios of virial to lumi- 



nous mass; c.f., Figurejl]). Thus, GMCs in dwarf galaxies 
are remarkably compatible with a Galactic Xqq inde- 
pendent of their metallicity, and any metallicity trends 
appear to be much weaker than the dispersion of the 
measurements . 

To further quantify this statement we have carried out 
a number of fits to the data. We do not consider in 
the fits NGC 2976 (because of its uncertain metallicity), 
and the Milky Way (which may suffer from a system- 
atic methodological offset). A least-squares bivariate fit 
log Xqo ~ a + b log{0/H) with simultaneous errors in 
both axes to the measurements for all remaining galaxies 
yields a slope b = —0.45 ± 0.30 and x^ ~ 21, with errors 
estimated using bootstrapping. Increasing all errors pro- 
portionally to obtain x^ « 1 yields b = —0.23 ± 0.25. 
Similar fits to the data for the galaxies with the largest 
number of GMC measurements (M 31, M 33, IC 10, 
LMC, and SMC) yield b = -0.46 ± 0.30 (x^ ~ 5) and 
b = -0.24 ± 0.26 (x^ « 1). These resufis, plus the con- 
siderations in the paragraph above, confirm that there is 
no measurable trend in the resolved Xqo with metallic- 
ity present in the data. As we discuss below, however. 
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this does not imply that there is no metallicity trend in 
the global ratio of CO-to-H2 in galaxies. 

What about the possibility of clouds being more mas- 
sive than the simple virial mass estimate? We already 
discussed the observation that several clouds in IC 10 
and the SMC are not quite compatible with the Galactic 
size-line width relation, in the direction of too small a 
line width for a given size. We will argue in ii5.4l that 
a possibility is that these clouds are transient structures 
not in equilibrium that will collapse in a few Myr due 
to the lack of turbulent support, or that they may be 
supported by magnetic fields in larger proportion than 
the typical Galactic GMC. In either case, such clouds 
could be removed from the sample since their line width 
may underestimate their true mass. Doing so slightly 
changes the averages for the SMC and IC 10, but does 
not significantly alter the results. 

Thus, a consistent and uniform analysis of the avail- 
able data for resolved GMCs shows no evidence for an 
increasing Xqo with decreasing metallicity. We note 
that this is unlike the scenario mentioned before, where 
similar gradients in radiation field and metallicity con- 
spire to keep Xco approximately constant in the disks 
of the Milky Way or M 33. Here we have galaxies (and 
particular GMC complexes) that simultaneously have a 
higher-than-Galactic radiation field and a lower-than- 
Galactic metallicity, yet a very s i milar Xgo for their 
resolved GMCs (e.g., [Israel et all Il996t iMadden etall 
[19971: IWilke eraI][20M^ 

We emphasize that these results correspond to the re- 
solved Xco- As we discussed in %.'2.'2[ FIR observa- 
tions show that in galaxies such as the SMC, CO-bright 
cores are likely embedd ed in extended envelopes of II2 
that do not emit in CO (jMadden et al.lll997l:lLerov et all 
l2007f) . The global molecular gas content of a galaxy rel- 
ative to its CO luminosity may well steeply scale with its 
dust-to-gas ratio or metallicity, as it is strongly su ggested 
by the available data o n a few of these objects (|Israell 
119971: iLerov et ahl 120071) . Nevertheless, it appears that 
the structures that we are able to identify as individ- 
ual GMCs by means of their CO emission have CO-to- 
H2 conversion factors (as well as Larson relations) that 
are approximately Galactic, independent of their nebular 
metallicities. 

We can understand the joint results from the CO and 
dust-continuum observations in the following terms: in 
a low-metallicity environment such as the SMC, CO- 
bright clouds are the innermost portions of consider- 
ably larger II2 structures mostl y devoid of CO emis sion. 
This scenari o was suggested bv iRubio et al.l (|1991h and 
iRubio et al.l ([l993b) for the SMC and is consistent with 
calculations of the effect of metallici ty on the placement 
of the photodiss ociation fronts (Malonev fc Black. 1988. : 
lElmegreenlll989t ). It also ha s been suggested as the sit- 
uatio n in outer gal axy disks ("Papadopou los. Thi. fc VitH 
I2OO2I but see also IWolfire e t al. 2008 for a recent de- 
termination of H2 formation rates in diffuse gas). It 
appears that (as de scribed for Milky Way GMCs by 
iHever fc Bruntj[200^ these inner portions approximately 
follow the Larson relations and, to CO observers capable 
of resolving them, they appear very similar to Galac- 
tic GMCs. Thus an approximately Galactic Xco factor 
correctly estimates the mass of the CO-bright core, but 
underestimates considerably the mass of II2 in the entire 



GMC. In this scenario the bulk of the mass has to be lo- 
cated in an envelope, surrounding the CO emission. This 
requires that the radiation field responsible for the pho- 
todissociation of CO be external to the cloud, and that 
II2 be effective at self-shielding even in an environment 
with low dust-to-gas ratio and consequently reduced II2 
formation rate. 

5.3. Brightness Temperatures 

Column (11) in Table[3]gives brightness temperature at 
the cloud peak, T^, for each GMC, a quantity set by the 
excitation temperature, optical depth, and filling factor 
of CO. Although not part of the Larson relations, this 
quantity is accessible in extragalactic GMCs and may 
yield insight into the physical state of clouds. 

For Galactic clouds the average brightness tempera- 
ture is ~ 4 K ([Solomon et al.lll987l ). Many of our data 
exhibit notably lower T^ than this. This might indicate 
either lower excitation temperatures or a different optical 
depth of CO, due to e.g. a lower cloudlet filling factor 
or other geometry well below our resolution. However, 
the simplest explanation is that we observe GMCs with 
a (spatially) large beam and that beam dillution, i.e. a 
low CO filling factor within the beam, lowers T^ in more 
distant systems. 

Figure [7] shows that beam dilution indeed explains 
many of the the low observed T^. We plot the mean 
and full range of T^ for each galaxy as a function of 
the spatial resolution of the data. We also show the ex- 
pected Tb for the case of a uniform brightness cloud 
with FWHM size 40 pc (dashed line). The dashed line 
is consistent with most of our data, including both Local 
Group galaxies (the LMC, IC 10, M 33, and M 31) and 
more distant dwarfs (NGC 2976, NGC 4214, NGC 4605). 
For the more distant dwarfs, this highlights that while we 
resolve large (^ 100 pc) structures in the CO, these are 
likely blends of several GMCs. 

A few extreme values of T^ cannot be immediately 
explained by beam dilution. NGC 3077 appears un- 
usually bright in CO when observed with both BIMA 
and OVRO, perhaps a sign the that ongoing starburst at 
the heart of this galaxy leads to higher excitation tem- 
peratures. NGC 4449 shows a similarly high T^. In 
NGC 185, NGC 205, NGC 1569, and the SMC (N83) Tb 
is low despite excellent spatial resolution. In these galax- 
ies, molecular clouds are physically smaller, have lower 
area filling fractions of CO emission, and/or have lower 
kinetic temperatures than in large galaxies. The first 
two might be caused by the low abundance of CO or by 
diminished shielding from dust as discussed in previous 
sections. A lower kinetic temperature remains a possi- 
bility, but is not suppor ted by existing studies at least i n 
the Magellanic Clouds (jBolatto. Israel, fc Martinll2005| ). 
Clouds in the SMC are indeed physically smaller than in 
most other galaxies, but that effect does not completely 
account for the observed differences (see dotted line in 
Figure [7|). This suggests that the CO emission has a 
lower area filling fraction: clouds in some of these galax- 
ies are porous structures, where CO arises only from high 
Av clumps. 

5.4. The departures in the Small Magellanic Cloud 

We discussed in Ji4.2l the peculiar situation of several of 
the least massive clouds in this study, mostly those be- 
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longing to the SMC, that systematically exhibit velocity 
dispersions that are too small for their sizes according 
to the Galactic size-line width relation. We noted that, 
despite their small sizes, these objects display exactly 
the opposite behavior from that observed in the outer 
Galaxy, where small clouds are confined by the external 
pressure (jHever et al.ll200l|) . 

Under the assumption of turbulence-supported clouds 
in virial equilibrium following a ct„ oc BP'^ relation, the 
fact that many SMC and a few IC 10 clouds lie a factor of 
two in Ur, below the size-line width relation implies that 
their surface densities are four times lower than Galac- 
tic clouds, yielding Sgmc ^ 45 Mq pc~^. Recall that, 
for Milky Way clouds, the observed surface density trans- 
lates into a visual extinction Ay ~ 7.5 through the cloud. 
In low-metallicity galaxies such as the SMC and IC 10, 
where the dust-to-gas ratio is lower than in the Milky 
Way by at least a factor similar to their heavy-element 
deficit, a reduction by a factor of 4 in surface density is 
compounded by another factor of ~ 4 in the dust-to-gas 
ratio implying that we would expect extinctions of order 
Ay ^ 0.5 through the cloud, or Ay ~ 0.2 — 0.3 at the 
cloud center. We do not expect bright ^^CO emission at 
such low extinction, where most CO molecules would be 
photodissociated. It is possible, however, that GMCs in 
the SMC have a considerably lower CO area filling frac- 
tion than their Galactic counterparts. Thus, although 
the average surface densities and corresponding extinc- 
tions are are too low to allow the formation of CO, it still 
exists in small, well shielded clumps within these struc- 
tures. We have discussed in the previous section that 
there is some evidence along these lines, since clouds in 
the SMC indeed have lower brightness temperatures sug- 
gestive of smaller CO area filling fractions. 



Alternatively, if the hypothesis of virialized clouds sup- 
ported chiefiy by turbulent motions is not valid, the SMC 
and IC 10 GMCs that fall under the size-line width re- 
lation show a deficit of turbulent kinetic energy with re- 
spect to similar size Milky Way clouds. Since Galac- 
tic clouds are supported against collapse by a combina- 
tion of turbulence and magnetic fields with turbulence 
lending the main support on the large scales (jMcKed 
119891 : iMcKee fc Ostrike"^ l2007t ) . this apparent deficit of 
turbulent energy would translate into rapid collapse and 
subsequent star formation in a free-fall timescale (i// ~ 
4.4 (n/lOO)^*^'^ Myr, where n is the mean volume den- 
sity of the cloud; lMcKedri999[ ) unless there is significant 
cloud support provided by a magnetic field that is pro- 
portionally larger (or better coupled to the cloud ma- 
terial, perhaps due to larger cloud ionization fractions) 
than in otherwise similar Galactic structures. 

To test this hypothesis it is necessary to obtain cloud 
mass determinations independent of virial assumptions. 
This is usually attained by employing the Xqq factor, 
modeling molecular line emission, or using measurements 
of the dust continuum and assuming a dust-to-gas ra- 
tio aiid_grain emissivity. The latter method was used 
bv lBot et al.l ()2007t ). who determined masses for several 
molecular clouds in the SW region of the SMC. These au- 
thors found the dust-derived masses for these clouds to 
be, on average, twice as large as their virial masses. Un- 
der the assumption of long-lived clouds, they attributed 
this fact to additional magnetic support by a B ~ 15 
/iG field. This field is similar to that present in smaller 
structures (a few parse cs in size with masses M < 10"^ 
Mq) in the M ilky Way (iBertoldi fc McKee 1992) . Along 
similar lines, iLerov et al.l ( 2007t) modeling of the far- 
infrared emission in the SMC to derive molecular gas 
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surface densities found a mean surface density of molec- 
ular gas Egmc — 180 ±30 M© pc~^ on 46 pc scales, very 
similar to that observed in Milky Way GMCs. These ob- 
servations suggest that surface densities in the SMC are 
higher than what would be implied by cloud line widths 
under the assumption of turbulence-supported virialized 
clouds. 

We argued in ij5.2.3l that the joint CO and dust con- 
tinuum data in the SMC is best understood in terms of 
extended CO-faint H2 envelopes surrounding CO-bright 
cores. In this interpretation, part (possibly a large part) 
of the mass or surface density obtained by dust contin- 
uum modeling is contributed by the extended H2 enve- 
lope, and the mass of this envelope may not be accurately 
reflected in the kinematics of the CO-bright core. Conse- 
quently, the high dust-derived surface densities and low 
observed CO velocity dispersions may not necessarily re- 
quire an enhanced magnetic support. We expect that 
further joint analysis of high spatial resolution FIR and 
CO observations will shed light on the existence of ex- 
tended II2 envelopes f|Lerov et al.ll2008[ ). 

6. SUMMARY AND CONCLUSIONS 

We present and discuss a comprehensive set of high 
spatial resolution observations of CO in low mass galax- 
ies, obtained by a combination of interferometer and 
single-dish instruments. Although the data are hetero- 
geneous, we analyze them in a consistent manner to ob- 
tain CMC sizes, velocity dispersions, and luminosities. 
Th e an alysis is per formed usin g the algorithm described 
by iRosolowskv fc Lerovl ()2006D . which does a good job 
at removing the biases due to dissimilar resolution and 
signal-to-noise. We compare this uniform dataset of re- 
solved extragalactic molecular cloud properties against 
those of GMCs in the three disk galaxies in the Local 
Group. T o do so, we analyz e the int erferometric maps 
of M 31 (|Rosolowskvll2007D . M 33 (jRosolowskv et all 
20031). and the samp le of Galactic GMCs discussed by 
Solomon et"al] ()1987[ ). 

The main result of this study is that, remarkably, there 
are only small differences between CO-bright GMCs in 
the Milky Way and GMCs in galaxies with metallicities 
as low as 0.2Zq subject to a variety of physical condi- 
tions. Our sample of extragalactic GMCs follows approx- 
imately the same size-line width, luminosity-size, and 
luminosity-line width relations as Galactic clouds. Such 
uniformity may be in part responsible for the observed 
invariance of the stellar Initial Mass Function in galaxies. 
In any case, this result underscores that the Galactic Lar- 
son relations provide a remarkably good description of 
CO-bright Giant Molecular Clouds independent of their 
environment, at least in the range of environments ex- 
plored by this study. 

Although the Larson relations are approximately 
Galactic there are some significant departures. GMCs 
in dwarf galaxies tend to be slightly larger than GMCs 
in the Milky Way, M 31, and M 33 for a given CO lu- 
minosity or velocity dispersion. The largest departures 
occur in the SMC, the galaxy with the lowest metallicity 
of the sample. A possible interpretation, viable for most 
of our objects, is that GMCs in small galaxies have on 
average a surface density of Sgmc ~ 85 Mq pc~^ rather 
than the canonical Sgmc — 170 M0 pc~^ observed in 
the Galaxy. In the case of the SMC, however, the im- 



plied surface densities, and consequently the central ex- 
tinction in the GMCs, would be too low to expect bright 
CO emission. We explore three possibilities: these clouds 
are transient structures not supported by turbulence and 
will collapse in a few Myr, or else they are supported by 
magnetic fields in a larger proportion than similar Galac- 
tic clouds, or maybe the kinematics of the CO cores do 
not reflect the presence of massive H2 envelopes. 

Analysis of the properties of our sample of extragalac- 
tic GMCs shows that they do not accord with simple 
predictions fro m the theory o f photoionization-regulated 
star formation (|McKedll9890 . In the framework of this 
theory, clouds in lower metallicity environments will have 
larger surface densities to attain similar extinction in 
their central regions, which can then decouple from the 
magnetic field and collapse. Such a trend would be ev- 
ident in the coefficient of the size-line width relation, as 
expressed by Equation [T3l As Fig. [5] illustrates, we see 
no evidence for such a trend in our data. We point out 
four possible caveats with these results: 1) it is possible 
that our observations do not probe the spatial scales on 
which these enhancements may occur. 2) Perhaps other 
parameters that enter in the theory (such as the cosmic 
ray ionization rate, for example) change from galaxy to 
galaxy in a manner that conspires to keep the coefficient 
of the size-line width relation approximately constant. 3) 
Maybe our identification of Sgr with metallicity is incor- 
rect and we sample a smaller range of conditions than 
we think we do. Or, 4) the CO kinematics do not trace 
the full potential of the cloud. In any case, we find a 
cr^ i?^"'^ product that is approximately constant or even 
decreasing for decreasing metallicity, contrary to expec- 
tations. 

Finally, we address the matter of the dependency of 
the CO-to-H2 conversion factor, Xco, on metallicity. We 
find that the extragalactic GMCs in our analysis agree 
very well with the Galactic luminosity-virial mass rela- 
tion. Consequently, in our study of resolved GMCs we 
find no measurable change of Xqo'- over a factor of 5 
in metallicity most of our galaxies are compatible with 
a Galactic Xco = 2 x 10^° cm~2(K km s^'^y^ within 
a factor of two, and there is no discernible trend with 
metallicity. We emphasize, however, that this measure- 
ment is relevant on the scales of the individual, CO- 
bright GMCs. Studies in the FIR and at millimeter- 
waves in low metallicity environments (particularly the 
SMC) suggest that these clouds are embedded in larger 
II2 envelopes that are not traced by CO emission but 
contain an appreciable mass of molecular gas, at least 
in the case of the SMC (jlsrael II1997I : lRubioetal.1 120041 : 
iLerov et"al]|2007l ). These FIR results are approximately 
consistent with virial Xco estimates obtai ned on larger 
scales than the ones discussed in this paper (iRubio et al.l 
ll993bHMi^uno et al.ll2001l : iBolatto et al.ll2003( l. We sug- 
gest that the FIR and CO observations can be simul- 
taneously understood in the following terms: in low 
metallicity gas bright CO emission is relegated to the 
density peaks. Observations that resolve those density 
peaks show properties that are similar to those of Milky 
Way GMCs — that is the result of this study, and it 
is supp orted by studies of th e structure of Milky Wa y 
GMCs (iHever fc BruntI [2OOI IRosolows kv eTall [200^. 
On larger scales, the cloud-to-cloud velocity dispersion 
as well as the dust opacity suggest that there is a large 
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mass of likely molecular material that is not accounted 
for by the CO intensity. Such regions would be CO-bright 
in objects of higher metallicity. 

Studying the resolved properties of extragalactic 
GMCs is a challenging undertaking with the current gen- 
eration of instruments, particularly for clouds in CO- 
faint dwarf galaxies. In the near future the deployment of 
the ALMA interferometer will make it possible to obtain 
more precise measurements on larger samples at farther 
distances. Studies such as the one presented here suffer 
from an important bias: we can only perform them on 
the brightest CO peaks. The detection of clouds as faint 
as those in the SMC in CO, for example, is currently 
impossible beyond the immediate neighborhood of the 
Local Group. These improvements in instrumentation 
will soon provide powerful tests of theories of star and 



molecular cloud formation beyond our own Galaxy. 
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